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ABSTRACT Cell-penetrating peptides (CPPs) such as HIV's trans-activating transcriptional activator 
(TAT) and polyarginine rapidly pass through the plasma membranes of mammalian cells by an unknown 
mechanism called transduction. They may be medically useful when fused to well-chosen chains of 
fewer than about 35 amino acids. I offer a simple model of transduction in which phosphatidylserines 
and CPPs effectively form two plates of a capacitor with a voltage sufficient to cause the formation of 
transient pores (electroporation). The model is consistent with experimental data on the transduction 
of oligoarginine into mouse C2C12 myoblasts and makes three testable predictions. 



I. CELL-PENETRATING PEPTIDES 

Polyarginine, TAT [HI2], Penetratin, and other short, 
positively charged peptides can penetrate the plasma 
membranes of live cells and can tow along with them 
cargos that greatly exceed the 600 Da restriction bar- 
rier. They are promising therapeutic tools when fused 
to well-chosen sequences of fewer than about 35 amino 
acids (3HI1]. 

TAT carries cargos across cell membranes with high ef- 
ficiency by at least two functionally distinct mechanisms 
according to whether the cargo is big or small [15] . Big 
cargos, such as proteins or quantum dots, enter via cave- 
olae endocytosis and macropinocytosis [T7] [TH] , and rel- 
atively few escape the cytoplasmic vesicles in which they 
then are trapped [16] . 

Small cargos, such as peptides of fewer than 30-40 
amino acids, enter both slowly by endocytosis and rapidly 
by transduction, an unknown mechanism that uses the 
membrane potential [T6lH9ti2T]. Peptides fused to TAT 
enter cells within seconds |22) . 

Se c. \H\ reviews some facts about plasma membranes. 
Sec. |III| reviews how the phospholipid bilayer prevents 
ions from crossing the plasma membrane. Sec. |IV| de- 
scribes a simple model of transduction in which CPPs 
on the outer leaflet and phosphatidylserines on the inner 
leaflet form a kind of capacitor with a voltage sufficient 
to favor the formation of transient pores (electropora- 
tion). Sec. |v| shows that the model is consistent with 
measurements made by Tiinnemann et al. [23] on the 
transduction into mouse myoblasts of oligoarginines and 
oligolysines carrying fluorophores of 400 Da. Sec. |VI| tclls 
how to test three predictions of the model. The paper 
ends with a short summary in Sec |VII| 

II. PLASMA MEMBRANES 

The plasma membrane of a mammalian cell is a lipid 
bilayer that is 4 or 5 nm thick. Of the four main phos- 
pholipids in it, three — phosphatidylethanolamine (PE), 
phosphatidylcholine (PC), and sphingomyelin (SM) — are 
neutral, and one, phosphatidylserine (PS), is negatively 



charged. In live cells, PE and PS are mostly in the cy- 
tosolic layer, and PC and SM in the outer layer j2H US] . 
Aminophospholipid translocase (flippase) moves PE and 
PS to the inner layer; floppase slowly moves all phospho- 
lipids to the outer layer [24] . 

Glycolipids make up about 5% of the lipid molecules 
of the outer layer of a mammalian plasma membrane. 
Their hydrocarbon tails normally are saturated. Instead 
of a modified phosphate group, they are decorated with 
galactose, glucose, GalNAc = N-acetylgalactosamine, 
and other sugars. The most complex glycolipids — 
the gangliosides — have negatively charged sialic-acid 
(NANA) groups. Incidentally, cholera toxin binds to and 
enters cells that display the Gmi ganglioside. [35] 

A living cell maintains an electrostatic potential of be- 
tween 20 and 120 mV across its plasma membrane. The 
electric field E within the membrane points into the cell 
and is huge, about 15 mV/nm or 1.5 x 10 7 V/m if the 
potential difference is 60 mV across a membrane of 4 
nm. Conventionally, one reports membrane potentials 
as the electric potential inside the cell minus that out- 
side, so that here AV = —60 mV. Near but outside 
the membrane, this electric field falls-off exponentially 
E(r) = E exp(— r/D() with the ratio of the distance r 
from the membrane to the Debye length Dg, which is of 
the order of a nanometer. The rapid entry of TAT fused 
to peptides is frustrated only by agents that destroy the 
electric field E [16] . which applies a force qE to a CPP 
of charge q. 

Most of the phospholipids of the outer leaflet of the 
plasma membrane are neutral PCs & SMs. They vastly 
outnumber the negatively charged gangliosides, which are 
a subset of the glycolipids, which themselves amount only 
to 5% of the outer layer. Imagine now that CPP-cargo 
molecules are in the extra-cellular environment. Many 
of them will be pinned down by the electric field E(r) 
just outside the membrane, their positively charged side- 
chains interacting with the negative phosphate groups of 
neutral dipolar PC & SM head groups. (Other CPP- 
cargo molecules will stick to negatively charged gan- 
gliosides and to glycosaminoglycans (GAGs) attached to 
transmembrane proteoglycans (PGs); these slowly will be 
endocytosed. Heparan-sulfate PGs are needed for TAT- 
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protein endocytosis 26J.) It is crucial that the dipolar PC 
& SM head groups are neutral and so do not cancel or re- 
duce the positive electric charge of a CPP-cargo molecule. 
The net positive charge of a CPP-cargo molecule will at- 
tract negatively charged PSs and cause them to diffuse 
toward the part of the inner leaflet directly below the 
CPP-cargo molecule. This is the starting point for the 
model described in the next two sections. 



III. THE PUZZLE 

The dielectric constant q s» 2 of the hydrocarbons of 
a lipid bilayer is much less than that of water e w 80. 
Thus, the difference AE w ^g in the electrostatic energy 
of an ion of charge q and effective radius a in the bilayer 
and in water [57J is 



AE„ 



1 



1 
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or 3.5 eV if the ion's charge is that of the proton and 
its radius is a = 1 A. This energy barrier is far larger 
than the 0.06 eV gained when a unit charge crosses a 60 
mV phospholipid bilayer. Thus, an ion will not cross a 
cell's plasma membrane unless a transporter or a channel 
facilitates (and regulates) its passage. 

In the present model of CPP transduction, the electro- 
statics of the CPP-cargo complex and the role of PSs on 
the inner leaflet play key roles. 

The electrostatics of a cationic polypeptide such as 
TAT or polyarginine are more complex than for an ion. 
I will model the CPP and its cargo in water as a sphere 
with its positive charges on its surface. The density of a 
protein of mass M kDa is estimated [55J to be 



p(M) = 0.8491 + 0.0873 e 



-M/13 



kDa/nm 3 . (2) 



A CPP-cargo complex would not be expected to fold as 
densely as a natural globular protein, and so for it the 
estimate p(M) is more of an upper bound. The radius 
r of a putative sphere consisting of M kDa of CPP and 
cargo then would be at least 
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For instance, a CPP of N arginines and a tiny fluorophorc 
cargo of 400 Da has a mass of M N = 0.15627V + 0.4 kDa, 
and so its radius would satisfy 



> 



3 M N 
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which gives r = 0.75 nm for N — 8 arginines. The lower 
bounds on the radii for N = 5-12 are listed in column 2 
of the Table ^. 

For larger cargos of N c = 50-100 amino acids of 130 



Da each, the lower bounds on the radii range from 1.25 to 
1.59 nm. (In what follows, iV c will represent the number 
of amino acids in the cargo or the mass of the cargo in 
Daltons divided by 130 Da.) Adding another 0.8 nm for 
the PC/SM head groups extends these lower bounds on 
the radii to 2.05-2.39 nm. The diameters of these spheres 
would approach or exceed the normal 4.4 nm of the lipid 
bilayer. The extra energy penalty due to the distortion 
of the plasma membrane may explain why CPPs can not 
transduce cargos of more than 50 amino acids [16] . 

If the CPP-cargo molecule were a charged conducting 
sphere of radius r and charge q, then its electrostatic 
energy in water would be 



E(N,N c ,q,w) 
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This term neglects the short-distance detail of the electric 
field near the q/e positive unit charges e of the CPP-cargo 
molecule. So a short-distance correction term 



E sdc (N,N c ,q,w) 



qe 



(6) 



proportional to q must be added to E(N, N c , q, w). The 
short distance a is an adjustable parameter, which should 
turn out to be several A since the term E s( i c is a correc- 
tion to be added to E(N,N c ,q,w) and not the entire 
electrostatic energy. 

As the CPP-cargo molecule enters the lipid bilayer, 
the phosphate groups of the PC and SM of the outer 
leaflet bind to the positively charged guanidinium and 
amine groups of the CPP. The positive charges of the 
phosphocholine groups of PC and SM are about d = 5 A 
from their phosphate groups [35] . The binding of PC and 
SM therefore approximately increases the effective radius 
of the charged sphere to r m ss r + d. The electrostatic 
energy of this complex in the hydrocarbon tails of the 
lipid bilayer then is 



E(N,N e ,q,£) 



8ne ei(r + d) 



apart from a correction factor 

E sdc (N,N c ,qJ) 
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(8) 



similar to Q. 

Apart from correction terms, the electrostatic energy 
penalty when the CPP-cargo molecule enters the lipid 
bilayer from water as a CPP-cargo-PC/SM complex is 
the difference 

AE° w>t {N,N c ,q) « E(N,N c ,q,£) - E(N, N c ,q,w)(9) 

r + d 6( 



&ire e e (r + d) 



Because the thickness t — 4.4 nm of the lipid bilayer is 
only a few times that of the CPP-cargo-PC/SM complex, 
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we also must include the Parsegian correction 



In 



2e„ 



(10) 



The sum of the water-to-lipid energy ^ and Parsegian's 
correction ( 10 1 is 



AE Wtt (N,N c ,q) = AE^ e {N,N c ,q) + AE P . (11) 

The energy AE W e(N, N c , q) is listed in column 3 of Ta- 
ble |l| for a CPP of N = 5-12 arginines towing a fluo- 
rophore cargo of 400 Da with d — 0.5 mn. 

The short-distance correction terms augment this 
penalty by 



AE sdc (N,N c ,q) 



E sdc (N, N c , q, I) - E sdc (N, N c , q, w) 

qe ( ^ e £ 
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and do not require Parsegian's correction because they 
are short-distance effects. This short-distance correction 
AE sdc is listed in column 4 of Table Q for a CPP of N = 
5-12 arginines and a representative value of a = 4.5 A 
for the short-distance parameter. 

The net electrostatic energy penalty when the CPP- 
cargo molecule enters the lipid bilayer from water as a 
CPP-cargo-PC/SM complex is then the sum of (|9j[l0j & 
12) 



AE„ 



AE 



AE P + AE, 



(13) 



A CPP of 8 arginines carrying a fluorophore of 400 Da 
has a radius r of 0.75 nm, and with a — 4.5 A, the change 
( 13 1 in its electrostatic energy on going from water to 



lipid is 



AE, 



,3,8e) « 16.9 eV. 



(14) 



This energy barrier is 35 times bigger than the energy 
0.48 eV that it gains by crossing a potential difference of 
60 mV. So how and why does it cross? 



IV. A CPP-PS CAPACITOR AND 
ELECTROPORATION 

My answer lies in the second part of the present 
model — the phosphatidylserines (PSs). They comprise 
some 8-18% of the inner leaflet by weight [3D]. The PSs 
diffuse laterally within that leaflet with a diffusion con- 
stant D « 10~ 8 cm 2 /sec [3T] and so within one second 
spread to an area of 12 /im 2 , which is a significant frac- 
tion of the surface area of a eukaryotic cell. Attracted 
by the positive charge of two or more CPPs, several PSs 
cluster in the inner leaflet below the CPPs forming with 
them a kind of capacitor. In what follows, I will focus on 
the PSs because they are constrained to lie in the two- 
dimensional inner leaflet while other anions drift in the 



TABLE I: The radius r of a CPP-cargo molecule of N 
arginines and a cargo of 400 Da, its change in electrostatic 
energy AE W> £ when transferred from water to hydrocarbon 
(111, and the short-distance correction AEsdc 1 12 1. Dis- 
tances are in nm and energies in electron- Volts. 



N 


r 




AE sdc 


5 


0.67 


4.61 


3.90 


6 


0.70 


6.39 


4.68 


7 


0.73 


8.41 


5.46 


8 


0.75 


10.64 


6.24 


9 


0.78 


13.06 


7.02 


10 


0.80 


15.68 


7.80 


11 


0.82 


18.48 


8.58 


12 


0.84 


21.44 


9.36 



three-dimensional cytosol. 

After oligoarginine CPPs are added to the extra- 
cellular fluid medium, the arginines of some of them will 
bind to the phosphate groups of the phosphatidylcholine 
(PC) and sphingomyelin (SM) phospholipids in the outer 
leaflet. The electrostatic potential of an isolated CPP 
may not be strong enough to form a pore because it is 
reduced at a point r in the lipid by the mean relative 
permittivity e — (e w + £f)/2 to V(r) = q/ (Aireoer) and 
further reduced by counterions and image charges. But 
if a second CPP or a few CPPs should land near the first 
CPP and bind to PCs or SMs, then transduction becomes 
possible in the present model. 

A CPP of N arginines can form an a-helix of length 
L fa 0.16 V nm (or a /3-strand of length 0.34 V nm). If 
two R N CPPs lie a distance d apart bound to phosphate 
groups of the outer leaflet, then the charge 2Ve of the two 
CPPs is spread over approximately L d, giving a surface- 
charge density of a — 2Ne/(dL). The electric field of 
the two oligoarginines would attract 2V PSs to the part 
of the inner leaflet below the two K N s forming a kind 
of capacitor. If the voltage drop across this capacitor 
exceeds the threshold for the formation of the pre-pores 
or pores of electroporation, then the CPP and its cargo 
may enter the cell. This threshold is between roughly 
—250 and 550 mV, depending upon the duration of the 
potential [32] • 

Due to its bidentate guanidinium group, arginine binds 
to phosphate groups better than lysine [2TJ[33]. Thus, 
oligoarginines would be more likely than oligolysines to 
bind to the phospholipids of the outer leaflet forming 
transiently stable capacitors with the PSs of the inner 
leaflet. 

CPPs with V = 12 arginines can bind to more than 
one phosphate group, and so they can bind to more than 
one phospholipid of the outer leaflet. The number of 
such binding sites drops with the number of Rs. In the 
"picket-fence model," transmembrane proteins bind to 
actin and restrict the motion of outer-leaflet proteins and 
lipids 34J; so I assume that oligoarginines with V > 8 
remain approximately stable for a few hundred ms. 
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To include the effects of entropy and of the geometry 
of the oligoarginines (R^s) and the phosphatidylserines 
(PSs), I wrote and ran a Monte Carlo program (included 
as supplementary material). In this code, I assumed that 
two R^s were parallel a-helices a distance d apart bound 
to phosphate groups of the outer leaflet. The PSs were 
allowed to move in two dimensions in the inner leaflet at 
sites r*fc for k = 1, . . . 2N. The 2N positive charges of the 
two R N s are taken to lie fixed at the points 



' jCPP 



= (±d/2, 0.16(J\T/2-j), t) (15) 



in which the ±-signs label the two R s, and j labels the 
Rs of the R N s ( < j < N). The electrostatic energy of 
the PSs then is approximately 



TABLE II: The voltage differences AV (mV) across the 
plasma membrane due to two oligoarginines on the outer 
leaflet separated by a distance d nm and a flock of 2N phos- 
phatidylserines fluttering below them on the inner leaflet. 



AT 


AT/ r l O ^ 

zA V a — I nm 


AT/ A Q ^ «~ 

IXV a = 6 nm 


Za V a — 4 nm 


5 


-342 


-229 


-173 


6 


-398 


-266 


-198 


7 


-449 


-300 


-223 


8 


-504 


-335 


-248 


9 


-553 


-369 


-273 


10 


-601 


-401 


-297 


11 


-646 


-433 


-319 


12 


-687 


-465 


-340 



N PS i-1 



E 



s EE 



qe 



4w \££*=^ r *- r *l 



N PS N 

EE- 

fe=i j= i ~ r Tcpp\ 



in which e is the mean permittivity 

e=^(e™+e*). (17) 

The Monte Carlo code used a simple Metropolis step in 
which the x-y coordinates of all the PSs were randomly 
varied, each by less than ±0.025 nm. The code accepted 
any change in the positions of the PSs that lowered their 
energy E as given by ( 16 ) and accepted any change that 
raised their energy by AE conditionally with probability 



p _ e -AE/(kT) 



(18) 



in which k is Boltzmann's constant and T is 37 Celsius. 
The simulation took place in a two-dimensional box with 
side L = U N/0.13 nm, and periodic boundary conditions 
were enforced to keep the number density of PSs constant 
at 13% within the box. 

Each simulation started from a random configuration 
of PSs on the inner leaflet and ran for at least 10 mil- 
lion Metropolis steps. After at least 2 million steps for 
thermalization, measurements of the positions of the PSs 
and of the electrostatic potential were taken every 1000 
steps, and the positions of the PSs were recorded every 
80,000 steps. The distnace d separating the two CPPs 
was taken to be d = 2, 3, and 4 nm. (The simulations 
for d = 3 used ten times as many steps.) 

The voltages were measured at the points (0, 0, t) and 
(0,0,-0.1) nm. The offset by 0.1 nm was to avoid the 
chance arrival of a PS at the point (0,0,0). The voltage 
differences AV = 1/(0,0,-0.1) - V(0,0,f) are displayed 
in the Table [TTJ At all three separations, CPPs with 9 
or more Rs produce voltages that exceed the threshold 
( ~ 250-~ 550 mV ) for the formation of pre-pores or 
pores [32 . CPPs with 10 or more Rs may exceed the 
threshold for irreversible electroporation. 

The mean value of the distance of the PSs from the 



point (0,0,0) varied between (r) = 11.8 nm for two R 5 s 
separated by 2 nm and (r) = 27.5 nm for two R 12 s sepa- 
rated by 4 nm. 

In the present model, a CPP-PS capacitor increases the 
membrane potential sufficiently to cross the threshold for 
electroporation, at least for R^s with enough arginines. 
In its use of an electric field and of the binding of the 
CPPs to the phosphate groups of the phospholipids of the 
outer leaflet, the model has something in common with 
the adaptive-translocation model of Rothbard, Jessop, 
and Wender [5T]; in its use of neutral dipolar PC & SM 
head groups it is somewhat similar to the work of Herce 
and Garcia [35] and of Tang, Waring, and Hong . The 
key distinctive feature of the present model is the role of 
PSs in forming with the CPPs a CPP-PS capacitor with 
a voltage high enough to cause electropration. 



V. COMPARISON WITH EXPERIMENT 

Tiinnemann et al. |23] measured the ability of the L- 
and D-isoforms of oligolysine and of oligoarginine to carry 
fluorophores of ^400 Da into live cells. They found that 
oligoarginines transduced the fluorophores much better 
than oligolysines, that more arginines meant faster trans- 
duction, with L-R9 and L-R10 doing better than their 
shorter counterparts, and that the D-isoforms worked 
better than the L-isoforms, probably because of their 
greater resistance to proteases. 

The present model is consistent with these experimen- 
tal facts and explains them as follows: The oligoarginines 
crossed cell membranes more easily than the oligolysines 
because they were better able to bind to the phosphate 
groups of the PCs and SMs in the outer leaflet; the 
oligolysines were not able to form a stable upper plate of 
the CPP-PS capacitor. CPPs with more arginines were 
transduced more rapidly because with more arginines 
they could bind to more PCs and SMs and because their 
higher charge led to a higher trans-membrane potential. 
The D-isoforms worked better than the L-isoforms be- 
cause the capacitor mechanism is insensitive to the chi- 
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rality of the amino acids and because proteases were less 
able to cut them. 

This consistency of the capacitor model lends it some 
plausibility. The simplicity of the model suggests that it 
may apply to other oligocationic CPPs. But evolution 
finds what works, not what fits neatly into a model, and 
so other CPPs with different cargos may enter different 
cells by different mechanisms. 

VI. HOW TO TEST THIS MODEL 

One way to test the model would be to compare the 
rates of polyarginine transduction in wild-type cells and 
in those that have little or no phosphatidylserine (PS) in 
their plasma membranes. If PS plays a key role as in the 
model of this paper, then the transduction of polyargi- 
nine fused to a cargo of less than 30 amino acids should 
be faster in the wild-type cells than in those without PS 
in their plasma membranes. Mammalian cell lines that 
arc deficient in the synthesis of phosphatidylserine do ex- 
ist [37MT] , but they appear to have normal levels of PS 
in their plasma membranes }41) — presumably due to a 
lower rate of PS degradation [32] ■ 

Another test would be to construct artificial asymmet- 
ric bilayers [43H45] with and without PS on the "cytoso- 
lic" side and to compare the rates of CPP-cargo trans- 
duction. If the present model is right, then the rate of 
transduction should be higher through membranes with 
PS on the cytosolic side than through membranes with 
no PS or with PS on both sides. 

If CPPs do enter cells via molecular electroporation, 
then it may be possible to observe the formation of the 
transient pores by detecting changes in the conductance 
of the membrane |32j . Such measurements would be a 
third test of the model and would let one determine both 
whether CPP-transduction is related to the presence of 
PS on the cytosolic side of the membrane and whether it 
proceeds via molecular electroporation. 



VII. SUMMARY 

Cell-penetrating peptides can carry into cells cargos 
with molecular weights of an order of magnitude greater 
than the nominal limit 500 of the "rule of 5" 46J. Sec- 
tion |IV| describes a model in which phosphatidylserines 
(PSs) play key roles in the transduction of CPP-cargo 
molecules. In this model, the positively charged CPP 
on the outer leaflet and the negatively charged PSs on 
the inner leaflet form a kind of capacitor with a trans- 
membrane potential sufficient for electroporation (-250 
to -550 mV) . The model is consistent with experimental 
data [23] on the transduction of the L- and D-isomers of 
oligoarginine and oligolysine CPPs into mouse myoblasts. 

The model predicts that mammalian cells that lack 
phosphatidylserine in their plasma membranes transduce 
oligocations poorly, that artificial asymmetric bilayers 
with PS on the cytosolic side transduce oligocations bet- 
ter than ones without PS, and that transduction is ac- 
companied by brief changes in the conductance of the 
membrane. 
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